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1 Introduction 

This paper deals with a problem of timely detection and monitoring of re-shaping environmental boundaries. 
This task is of interest for tracking of oil or chemical spills m, exploration of radioactively contaminated areas, 
and for many other missions, such as tracking forest fires [S], or contaminant clouds m, or harmful algae 
blooms Eg, exploration of sea temperature and salinity or hazardous weather conditions like tropical storms, 
etc. In fact, such missions are devoted to localization of a spatially distributed phenomenon whose boundary 
can typically be defined as the level set (isoline) where a certain scalar field (e.g., the radiation level or the 
concentration of a pollutant) assumes a critical value. In many such missions, sensors observe the phenomenon 
only in a point-wise fashion at their locations. 

Networks of static sensors typically require high both deployment density and computational/communication 
loads to provide a good accuracy of observation [27]. More effective use of sensors is achieved in mobile networks, 
where each sensor can explore many locations and the sensors may be driven to the explored boundary, thus 
concentrating the overall sensing capacity of the network on the structure of the main interest. 

Autonomous unmanned mobile robots and robotic networks have been greatly used in the recent past for a 
variety of such missions in hazardous and complex environments not only to avoid or mitigate risks to human 
life and health but also due to their lightweights, inexpensive components, and low power consumptions, see 
e.g. [Tlf7ll8inTlll21ll41ll8l29B36l37l and references therein. These robots are often subjected to limitations on com¬ 
munication and so should be equipped with navigation and control systems that enable them to autonomously 
operate either constantly or for extended periods of time and distance. 

Recently, designs of such systems have gained much interest in control community. Many works in this area 
assume access to the field gradient or even higher derivatives; see, e.g., EunmEg and literature therein. 
For example, they present gradient-based networked contour estimation .IZlEIlGin.i centralized control laws 
stemming from the ’snake’ algorithms in image segmentation EH, potential-based approach m , estimation 
of the gradient and Hessian of a noisy field for driving the center of a rigid formation of sensors along a 
level curve [38] . However, spatial derivatives are often not directly measured. Meanwhile gradient estimation 
needs tight-flocking concentration of the sensors near the examined point, contrary to efficient utilization of the 
network, which ideally means its uniform distribution over the explored boundary. Data exchange needed for 
this estimation may be critically degraded by limitations on communication. Finally, sensor-based derivative 
estimators are prone to noise amplification, and their implementation is an intricate problem in practical setting 
HUSU3H- This carries a threat of performance degradation, puts strong extra burden on controller tuning, and 
may drastically increase the overall computational load [2j. 

A single mobile sensor with a point-wise access to the field value only is the main target for gradient-free 
approaches; see e.g., QUBUZHS and literature therein. Switches between two steering angles, which are carried 
out whenever the field reading crosses the desired value, were advocated in [161139] : a similar approach with more 
alternative angles was applied to an underwater vehicle in [5]. Such methods result in jagged trajectories and 
rely in effect on systematic sideways maneuvers to collect enough data, which gives rise to concerns about waste 
of resources. A method to control an aerial vehicle based on segmentation of the infrared images of the forest 
fire was presented in [8]. These works are based, more or less, on heuristics and provide no rigorous justification 
of the control laws and guarantees of success. In [B], a linear PD controller fed by the field value was proposed 
for steering a Dubins-car like vehicle with unlimited control range along a level curve of a radial harmonic field, 
and was supplied with a local convergence result. A sliding mode controller for tracking environmental level 
sets without gradient estimation was offered in EH- 

On the whole, previous research was confined to only steady fields for which the environmental boundary 
does not change over time. On the contrary, such boundaries are almost never steady in the real world, due 
to dispersion, advection, transport by wind, drift by water current, etc. Meanwhile the theory of tracking 
dynamic level sets lies in the uncharted territory. As a particular case, this topic includes navigation of a 
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mobile robot towards an unknowingly maneuvering target and further escorting it with a desired margirQ on 
the basis of a single measurement that decays away from the target, like the strength of the infrared, acoustic, or 
electromagnetic signal, or minus the distance to the target. Such navigation is of interest in many areas 
it carries a potential to reduce the hardware complexity and cost and to improve target pursuit reliability. To 
the best of our knowledge, rigorous analysis of such a navigation law was offered in [25] for only a very special 
case of the field — the distance to a moving Dubins-like car. However the results of [55] are not applicable to 
more general dynamic fields. 

This paper is aimed at filling the above gaps concerned with dynamic fields, while disembarrassing boundary 
tracking control from the intricacies related to gradient estimation and sideways fluctuations. To this end, we 
consider an under-actuated non-holonomic Dubins-car type mobile robot. It travels with a constant speed 
over planar curves of bounded curvatures and is controlled by the upper limited angular velocity. The tracked 
boundary is given by a generic dynamic scalar field. The robot has access to only the field value at the current 
location and the rate at which this reading evolves over time via, e.g., numerical differentiation. 

The proposed navigation law is non-demanding with respect to computation and motion. It develops some 
ideas set forth in [5B|; e.g., gradient estimates and systematic exploration maneuvers are not employed. Whereas 
the results of [Ml are not applicable to unsteady fields, the objective of this paper is to show that the ideas 
from [26] basically remain viable for generic dynamic fields. 

The extended introduction and discussion of the proposed control law are given in the paper submitted 
by the authors to the IFAC journal Automatica. This text basically contains the proofs of the technical facts 
underlying justification of the convergence and performance of the proposed algorithm in that paper, as well 
as illustrations of the main theoretical results, which were not included into that paper due to the length 
limitations. To make the current text logically consistent, we reproduce the problem statement and notations. 


2 Problem Setup and the Control Algorithm 


A mobile robot travels in a plane with a constant speed v and is controlled by the time-varying angular velocity 
u limited by a given constant u. There is an unknown and dynamic scalar field D(t,r ), i.e., a quantity D £ R 
that depends on time t and point r £ R 2 in the plane. Here r := (x,y) T stands for the pair of the absolute 
Cartesian coordinates x,y in R 2 . It is required to drive the robot to the spatial isoline D(t,r) = do where the 
field assumes a given value do and to subsequently drive the robot in a close proximity of this isoline so that the 
robot circulates along it. The robot has access to the field value d(f) := D(t , x, y) at the robot’s current location 
x = x(t),y = y(t) and to the rate d(t) at which this measurement evolves over time t, via, e.g., numerical 
differentiation of d(t). However, neither the partial derivative D ' t , nor D ' x , nor D' y is accessible. 

The kinematic model of the robot is as follows: 


x = vcos9, 9 = u , x(0) = x in . _ „ 

y = v sin 9, \u\<u, 2 /( 0 ) = y in ’ _ 1 


( 2 . 1 ) 


where 9 is the orientation angle of the robot. These equations are classically used to describe, for example, 
wheeled robots and missiles (see e.g. [22 [1231f5Tj and the literature therein). 

It is required to design a controller such that D[t,x{t),y[t)] —>■ do as t —)■ oo. 

In this paper, we examine the following navigation law: 

u{t) = — sgn {d(t) + x[d(t) - d 0 ]}u, (2.2) 


where sgn a is the sign of a (sgnO := 0 ) and x(') is a linear function with saturation: 


/ x IIP if \P\ < $ c , 0 ox 

x(p) ■■= \ XX ,, . . y-'yd. (2.3) 

I sgn (jp) \x otherwise 

Here 7 > 0 and S > 0 are design parameters. 

For the discontinuous controller (12.21) . the desired dynamics [Ml is given by d(t) = ~x[d(t) — do]. Since it is 
unrealistic to desire large d, saturation is a reasonable option. 


3 Quantities characterizing dynamic fields 

To judge the feasibility of the control objective and to tune the controller, we need some characteristics of the 
unsteady field D(-) or their estimates. Now we introduce the relevant quantities, along with employed notations. 

1 Which should be given by the desired field value. 
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• r a point on the plane; 

• (•; •} — the standard inner product in the plane; 

• !Rg = ^ — the matrix of counter-clockwise rotation through angle /3; 

• ^ = (/ib Sy) — the spatial gradient; 

• D" — the spatial Hessian; 

• /(i, d*) := {r : D(t,r) = d*} — the spatial isoline; 

• [T, N] = [T{t,r),N{t,r )}—the Frenet frame of the spatial isolinc I[t, d*] with d» := D{t, r), i.e., N(t, r) = 
||VDp’r)|| an< ^ un ^ t an g en f vector T(t,r) is oriented so that when traveling on I[t,d *] one has the 
domain of greater values := {r' : D(t,r') > d*} to the left; 

• x — the signed curvature of the spatial isolineH 

• ' \At\t, r) — the nearest (to r) point of intersection between the ordinate axis of the Frenet frame 
[T(t, r), N(t, r)] and the displaced isoline I[t + At, d*(f + At)], where the given smooth function d*(-) is 
such that d*(t) = D{t,r); see Fig. [T] 

• p d ^'\At\t, r) — the ordinate of {At\t, r); 

• A d *(')(i,r) — the front velocity of the isoline: A d *('\t,r) := limAt-s-o — () At t ^’ r ^ ') ^ d*(-) = const, the 
index is dropped in the last three notations; 


a(t,r) — the front acceleration of the spatial isoline d[t,d*],d* := D(t,r) at time t at the point r : 

A [t + At, r+(Ai)] - A[t, r] 


a(t,r) := lim 
v ' At->0 


At 


(3.1) 


where r+(At) := r + (At|t, r); 

• A(p(At\t, r) — the angular displacement of T[t + At, r + (Af)] with respect to T[t, r]; see Fig. [Q 

• u/(i,r)— the angular velocity of rotation of the isoline /[t,d*],d* := D(t, r), i.e., u(t, r) := limAt-»o ^ ; 

• q(Ad\t,r) — the ordinate of the nearest (to r) intersection between the N- axis of the Frenet frame 
[T(t,r),N(t,r)] and the isoline I(t\t, d* + Ad); 

• p(t,r) — the density of isolines at time t at point r: 


p(t,r) := lim 


Ad 


Ad—>o q(Ad\t, r) 




• v p (t,r) — the evolutional (proportional) growth rate of the above density at time t at point r: 


:= lim 


Ai t —^0 


p[t + At,r + (At)\- p(t,r) ' 
p(t,r)At 


(3.2) 


T p (t, r) — the tangential (proportional) growth rate of the isolines density at time t at point r: 


T p (t,r) := lim 


p(t,r +TAs) — p(t, r ) 


As— >o p[t, r)As 


(3.3) 


n p (t,r) — the normal (proportional) growth rate of the isolines density at time t at point r: 


,{t,r) := lim 


As—f0 


p(t, r + NAs ) - p{t, r) 
p{t,r) As 


(3.4) 


Now we show that under minor technical assumptions, these quantities are well-defined, and offer relation¬ 
ships that explicitly link them to D(-). 

2 This is positive and negative on the convexities and concavities, respectively, of the boundary of G^ . 

3 It assesses the “number” of isolines within the unit distance from I(t, d*), where the “number” is evaluated by the discrepancy 
in the values of D(-) observed on these isolines. 
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Lemma 3.1. //£>(•,•) is twice continuously differentiable in a vicinity of(t,r) andS7D(t,r) ^ 0, the quantities 
from Sec. [3] are well-defined and the following relations hold at ( t,r ) 


A = — * 


\ d *U = X + 


d* 


P=I|VD||, 


livin’.livin’ 

r + (dt) := r + (dt\t 1 r) = r + XNdt + o(dt ), 

(X7D' t + XD"N ; N) (VD! + A D"N: T) 

v 0 = -=- 


I V£>| 


'tt t A (VD' t ;N) A 
a = -- Xv 


D’! 




(D"T- T) ( D"N:T) 

x =- ,t 0 = , n„ = 


\X7D\\ 

pi 

(D"N-N) 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 


l|V£>|| ’ P ||VD|| 5 p ||VD|| 

Proof. The first claim and (13.51) . (13.61) follow from the implicit function theorem [IJj]. With (13.51) . (13.61) in mind, 
S7D [t + dt , r+(dt)\ = X7D[t + dt,r + XNdt + o(dt)\ = VD + [VD' t + XD"N]dt + o{dt), 


p(t + At,r+{dt)) H ||VU [t + dt, r + (df)]|| = || V£> + [VDj + XD"N]dt + o{dt)\\ 


= IIVDII 


(VD; VO; + A D"N) 

WM 


dt + o(dt) = || VD|| + (IV; VD' t + XD"N) dt + o(dt), 


(3.10) 


which gives the first formula in (13.71) . Furthermore, 


. , , , ,, VD\t + dt,r + (dt) 1 v 

N[t + dt , r + (dt)\ = -—^7 = N + o(dt) + 


\\X7D[t + dt,r + (dt)]\\ 


VDI + X D"N VD 


l|VD|| 


livin' 


(VD; S7D' t + A D"N) 


dt 


= N + 


1 [VD' t + XD"N - N (N: VD' t + XD"N)) dt + o(dt) = N + Tdt + o(dt), 


||V£>|| L - t ■ -- -■ -• v. ■ ^ ” 1 DVD,, 

where (a) holds since W = {W,T)T + (W,N)N \/W € K 2 . In the Frenet frame (T,N), we have 


N[t+dt,r + (dt)} = ( 


= N 


— cosO 

— sinO 


udt + o(dt) = N — uiTdt + o(dt). 


By equating the coefficients prefacing Tdt in the last two expressions, we get the second formula in (£3. 
Furthermore, 


\(t + it,r+(dt))^ - + <*,••+(<«)] E3J 


A- 


\\VD(t + dt,r + (dt))\\ 

D'{ t dt + {VD' t ;.r + (d.t)-r) , . (N;ND' t + XD"N) 


|VD|| 


D' 


IVDII 


dt + o(dt) 


m X _ D 


« - A ™ N) dt - X {N ^ D l + J D " N) dt + o(dt) 


l|VD|| 


|VD|| 


(b) a _ BlL 


-\(VD' t -,N) 

l|VD|| 


dt — Xvpdt + o(dt ) => (13.81) . 


where (b) follows from the first formula in (13.71) . The first equation in (13.91) is well known, the second one follows 
from the transformation 

(D"T:S7D) 

p{t,r + Tds) = ||V.D[t, r +Tds]|| = p(t,r) H- ||V7 ' r>|| — -ds + o(ds) = p(t,r) + ( D"T\N) ds + o(ds), 


l|VL>|| 


the third equation in (13.91) is established likewise. 


□ 
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4 Assumptions and the Main Theoretical Result 

To simplify the matters, we suppose that the operational zone of the robot is characterized by the extreme 
values d- < d+ taken by the field in it 

M := {( t,r) : d- < D{t,r) < d + } (4.1) 

and contains the required isoline d_ < do < d+. 

An extended discussion and justification of the following assumptions and results are presented in the main 
text supported by this paper. 

Proposition 4.1. Suppose that the moving robot remains on the isoline D[t,r(t)] = do, in a vicinity of which, 
the field is twice continuously differentiable and VD(-, •) 0. Then at any time the front speed of the isoline at 

the robot’s location does not exceed the speed v of the robot 

|A[f,r(t)]| < v, (4.2) 

the robot’s velocity v = ve,e = (cosd, sin$) T has the form 

v = XN ± vtT, where vt '■= \A 2 — A 2 , (4.3) 


and the following inequality is true 

|± 2 w + avf 1 + xvt | < u. ( 4 - 4 ) 

In ±, the sign + is taken if the robot travels along the isoline so that the domain {r : D(t,r) > do} is to the 
left, and — is taken otherwise. 

Assumption 4.1. The field is twice continuously differentiable in the domain (14.111 and there exist 

constants A\ > 0 and A u > 0 such that the following inequalities hold everywhere in AA: 


|A| < v - A a , 


(4.5) 


±2w + / . 9 + v 2 - A 2 


<u - A u , 


(4.6) 


where (1X61) is true with the both signs in ±. 

Since the robot’s path is smooth, it is reasonable to exclude isoline singularities. A common guarantee of 
their nonoccurrence is that the field has no critical points 132 - 

Assumption 4.2. At any time, the zone (14.11) does not contain critical points r of the field, i.e., such that X7D = 
0. Moreover, this property does not degrade: there is b p > 0 such that p(t,r) = ||VZl(t,r)|| > bf 1 V(t,r) £ AA. 

The next assumption is typically fulfilled in the real world, where physical quantities take bounded values. 

Assumption 4.3. There exist b\, b T , b n , bb v , b a € R such that the following inequalities hold in the set (HUD-- 


|A| < b\, 1 ' 7 -pi < b T , \n p \ < b n , |>r| < b\v p \ < b v , |a| < b a . (4.7) 

Under the control law m , the robot moves with u = ±u during an initial time interval. This motion is 
over the initial circle C±, which is the respective path from the initial state given in m- The last assumption 
requires that the encircled closed discs (also called initial ) lie in the operational zone (14.11) , and the maximal 

turning rate of the robot exceeds the average angular velocity of the field gradient, at least on some initial time 
interval. 

Assumption 4.4. There exists a natural k such that during the time interval [0,Tfc] ,7*, := (a) the gradient 

V D(t,ri n ),ri n := ( Xi n ,yi n ) T rotates through an angle that does not exceed 2n(k — 1 ) and (b) the both initial 
discs lie in the domain (14.11) . i.e., [0,Tfe] x D± C AA. 

Theorem 4.1. Let Assumptions 14.11 44.41 hold and the robot be driven by the navigation law ( 12 . 21 ) whose 
parameters 7 and p = yd satisfy the following inequalities, where /r* := b p p and a(pff) := yj A 2 — 2vp* — ffl: 


0 < /x* < 



A u > 



bx + 2 b v + 7 + 6 n /x* 


+ 


b a \ 


p*. 


(4.8) 


Then the robot achieves the control objective d(t) 


t—y OO , 

- > d 


0 ; 


always remaining in the domain m- 
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By sacrificing conservatism, the requirements (14.81) to the controller parameters can be transformed into an 
explicit form, which decrypts how small they may be. For example, let us pick £ £ (0,1) and confine the choice 
of (^*, 7 ) to the “semi-strip” 0 < /z* < ^Jv 2 + — v,j > 0 to ensure the first two inequalities in (14.81) . Then 

cr(/x*) > A\y/1 — ( and so (14.81) holds whenever 



bn + 2 b v + 7 + b n p^ 

VR 


A|(i-C ) 3 / 2 


/r* ^ A u . 


So recommended parameters lie in the portion of the above semi-strip that is cut off by the upper piece of the 
hyperbola explicitly described by the last relation provided that < is replaced by = in it. 

Remark 4.1. The last claim of Theorem 14.II implies that the behavior of the field outside the zone (14.11) does 
not matter. Moreover, Theorem 14.11 remains true even if the field is defined not everywhere outside (ED, which 
is of interest for some theoretical fields like c||r|| _1 . 

The fact that (14.81) holds for all small enough S provides a guideline for experimental tuning of the controller. 
To analytically tune it, estimates of the field parameters concerned in (14.81) should be known. Various types of 
knowledge about the field give rise to various specifications of (14.81) . The next section offers relevant samples. 


5 Some Particular Scenarios 

5.1 Dynamic Radial Field with Time-Invariant Profile 

Let it be known that the field is radial D(t,r) = c/(||r — r°(f)||) with time-invariant profile / : (0, 00 ) —> R. 
Then spatial isolines are circles with radii determined by the associated field levels d*. This opens the door for 
analytical computation of many field parameters from Section[3]and underlies transformation of (14.81) to a more 
explicit form. This may be used to acquire fully analytical conditions of isoline tracking for various special field 
profiles, which will be illustrated in the next two subsections. 

Let the profile / be twice continuously differentiable and decay f'(z) < 0; we initially assume that / is 
known. Conversely both field “intensity” c > 0 and its moving center r°(-) are unknown but obey some known 
bounds: 

0<R<||r in -r°(0)|| <R ||r°(t)|| < v 0 , ||r°(f)|| < a 0 , 0 < c_ < c < c+. (5.1) 

The objective is to advance to the held isoline I with the given level do and to subsequently track I. 

The control objective is well posed in the face of uncertainties only if the desired held value do lies in the held 
range for any c £ [c_,c+], or equivalently do/c± £ Im/. The following specihcation of Proposition 14.II shows 
that irrespective of the control law, the robot should be maneuverable enough to achieve the control objective, 
and details the required level of maneuvearbility. 

Proposition 5.1. Suppose that the robot is able to remain on the dynamic isoline D[t,r(t)] = do if the field 
uncertainties satisfy ED- Then the robot’s speed is no less than the maximal feasible speed of the field center 
v > vo, and the following inequality holds for any v n £ [0,uo], where R- := min{/~ 1 (do/c+); / -1 (do/c_)} and 
/ _1 (-) is the inverse function, 


a 0 


+ 


( y/v 2 -Vl+ V V 0 ~ V n) 


< U. 


\A 2 - Vl ’ R-yJv 2 - V 2 

Proof: An elementary calculus exercise, along with Lemma Id. 11 shows that 

W = -|jRj|- r=-* 5 JV, \ = {r°-,N), 

1 f" 

p = c \f'l *= || r _ r o|| > T p = 0 ’ n p = \fy v p = 0> 


(5.2) 


<r°;T) 


, a = (r°:N) + 


<*°;*r 


. (5.3) 


By Proposition 14.11 (14.21) and (14.41) hold. Putting v n := ( r°;N ) ,R := ||r — r°|| and noting that = 
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± \J ||f ,0 || 2 — v 2 , we see that now m •O- \v n \ < v and ED takes the form 


— u < —2 


|r°|| 2 -u 2 


(r°-N) + ^ 


R 


R 


'Jv ' 1 - u 2 


+ -bVv 2 - 


R 


|r°|| 2 -u 2 


R 


(r°;N} + Jl—1 

\/w 2 - vl 


+ - v 2 < u. 


This holds for any r°,r° satisfying the bounds from Ell). Minimizing and maximizing over feasible r° shows 
that 


- u < -2 


|f°|| 2 -u 2 


n —CiQ H~ 


R 


\/ v 2 -' 


+4 V v2 - u 

K 


\r°\\ 2 ~v 2 ap + ^M 1 


R 


\/v 2 - ■ 


+ V v2 ~ V n < ' 


R 


Putting r° := Vo N into \v n \ < v assures that vq < v. Minimizing and maximizing over feasible r° yields that 


Wv'o -vl , Qo ± v 2 +vp- 2 v n < - 
R y/v 2 - v 2 Ry/v 2 - v 2 

Since the expression following ± is nonnegative, the sign — can be dropped in ±. It remains to note that 
as c ranges over the feasible interval from ED. the parameter R runs over the interval with the end-points 
f~ 1 (do/c-) and f~ 1 (do/c + ) due to the isoline equation do = cf(R). □ 

The following corollary of Theorem 14.11 shows that slight enhancement of the above necessary conditions is 
enough for the controller ( 12 . 21 ) to succeed. Specifically, we in fact sharpen inequality (15.21) from < to < and 
extend it from the isoline on the transient by reducing f?_ in (15.21) . 

Proposition 5.2. Let the robot be faster than the field center v > vq and initially far enough from it: 


R in > (2v + Airvo)/u. 


(5.4) 


Suppose also that for some A u > 0 . ED holds for all v n G [0, «o] with u replaced by u — A u and R- given by 

R± := ± max{±i?(|_; ±/ _1 (d 0 /c_) ;±/ _1 (do/c+)} , 

R± := Rf n ± (2v + 4 ttv 0 )/li¬ 


nen there exist parameters 7 ,5 such that whenever ED holds, the controller ED brings the robot to the 
desired isoline: d(t) —> do as t — > 00 . Specifically, this is true if 


0 < /x* := 


7(5 


c- min \f(R)\ 
Re[R-,R+] 


< V y2 + i v ~ v o) 2 - V, 


where 


cr(p*) ■■= y/{v- vo) 2 - 2vp* - p 2 , 


■^7 + 7 + b n p* ao + 


<x{p* 


cr(p *) 3 


P* < A u , (5.5) 


\rm 

Re[R-fR + } \f'{R)\ ' 


b n = max 


Inequalities (15.511 are feasible and can be satisfied by taking small enough p := 71 5. Full knowledge of the 
field profile /(•) is not in fact required: it suffices to know estimates R°_ < R < R\ of the radius R of the 
desired isoline and replace i?_ and R + by their lower I?_ := minlJ?*).; R'_} and upper R + := maxjl?^; R' + } 
bounds, respectively. Furthermore, it suffices to know a lower minR e p{_,/?+] | f'{R)\ > 6 - > 0 and upper 
max^j g [ fl _ fl+ ] |/"(f?)| < b+ estimates and take p* = pcZ 1 bZ 1 and b n := b+/b- in Proposition l5.2l 

Proof of Proposition ED Let d'± := cf{R'^),d- := min{<f_; do}, := maxjd^; do}- Then due to (14.111 . 


M C {(f, r) : R- < ||r - r°(t)\\ < R + }. 


(5.6) 


Thanks to (15.311 . Assumptions 14.21 and l4. 31 arc fulfilled with b p = cl 1 ! min Re [ R _ /{ + ] |/ , (f?)|} 1 . As t ranges over 
[0, T 2 ], T 2 = Air/u, the field center r°(t) remains in the disc Dq of the radius Airvo/u centered at r°(0). This disk 
does not contain r in due to (15.41) and the first inequality in (15.11) . Meanwhile, the gradient VH(f,r in ) rotates 
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through an angle that does not exceed 7 r thanks to the first equation in (15.31) . and so the first requirement of 
Assumption 14.41 is satisfied with k = 2. Since D± C {r : ||r — r in || < 2 v/u}, for t £ [0,T 2 ],r £ D±, we have 
owing to m, 

||r - r°(f)|| > ||r in - r°(0)|| - \\r - r in || - || r°(t) - r°(0)|| > R^ n - 2v/u- Attv 0 /u = R'_, 

and similarly ||i— r°(t)|| < R' + . So the second claim of Assumption 14.41 with k = 2 holds due to (15.61) . 

Thanks to (15.31) . inequality (Id.51) is true with A a := v — vq, whereas (14.61) means that u — A u is no less than 


± 2 - 


<r°;T) (r°:N)+<0£ l 

nn I /-— 


r — r u 


yJv 2 ~(r°;N) 2 H r - r °H 


\] v 2 — (f°;2V)" 


Retracing the arguments from the proof of Proposition 15.11 shows that for (i, r) £ Ad, the last expression does 
not exceed that in (15.21) . Hence (Id.61) holds and Assumption 14.11 is fulfilled. 

Because of (15.31) and (15.61) . now the estimates (14.71) hold with 


b\ = vq, b T = 0 , b n = 


I/"WI 
* e “ aX fl + ] \f{R)\ ’ 


1 v 2 

b^ — ~r, : b v — 0 , b a — cio + ~r t ■ 
H— li¬ 


lt remains to note that this transforms (14.81) into (15.51) and apply Theorem 14.II □ 

The following remark shows that by sacrificing conservatism, the key condition (S3 can be simplified. 


Remark 5.1. Inequality (15.21) is true whenever 


ao . (v + vq ) 2 

— . 4-—- < u. 

y/v 2 — Vq vR- 

This holds since the addends in the left-hand side are upper bounds for the respective addends in (15.21) . 
Now we discuss more particular scenarios, where the tracking conditions can be further specified. 


N 



Figure 2: Escorting a maneuvering target. 


5.2 Approaching and escorting an unknowingly maneuvering target based on 
range-only measurements 

A pursuer robot should approach an unknowingly maneuvering target r°(t) using only measurements of the 
relative distance between them and should subsequently follow the target at the pre-specified range Rq, while 
always maintaining the constant speed v: see Fig. [2] The bounds from (15.11) are still valid; nothing else is known 
about the target. 

This scenario falls under the framework of subsect. 15.11 D[t, r} = — ||i— r°(t)\\,c = c± = 1, f(z) = —z, and 
do = —Rq- By Proposition 15.11 the pursuer is capable of carrying out the mission only if its speed is no less 
than that of the target v > vq and (15.21) holds with i?_ := Rq. Let these minimum requirements be fulfilled in 
the enhanced form given by Proposition [521 where i?_ := min{i? 0 : ftj/ — (2u + 47ruo)/?l}, and let (15.41) be true. 
By Proposition 15.21 the stated objective is ensured by the controller (12.21) if its parameters satisfy (15.51) . where 
/j* = 7(5 and b n = 0. 

In [25j . such problem was solved in the case where more information about the target is available: the target 
is a kinematically constrained Dubins-like car m with a known range u of angular velocities that moves at a 
constant speed obeying a known upper bound. 

















5.3 Detection and tracking of the boundary of a contaminated area transported 
by advection 

A planar liquid environment is chemically contaminated; see Fig. [3] The concentration D of the pollutant 
varies over time and plane, the critically contaminated area A cc is that where D exceeds a certain level do > 0. 
A robot should advance to the boundary of A cc and to subsequently track it, thus displaying A cc . Only the 
concentration at the robot’s current location can be measured. 

|| r — t»q || ^ 

The initial distribution of the concentration is Gaussian D(0,r) = ce ^ . Its center ro, variation cr 2 , 

and maximum value c are uncertain, but obey known bounds: 0 < cr_ < a < cr+,c > c_ > do ,0 < f?;“ < 
||r in — r 0 || < .fid. The liquid moves with a constant velocity V obeying a known bound ||Vj| < P*. Dispersion 
of the pollutant is negligible and its transport is basically by advection. 

Since D' t = — (V; VD) [2U] and so D(t, r) = D( 0, r—Vt ), this scenario falls into the framework of subsect. 1-5. II 
with r°(f) := r 0 +Vt,f(z) = e~ z d 2a \ao = 0,i>o = P*. By Proposition l5.il the robot is capable of carrying 
out the mission only if it is faster than the flow v > P* and (15.21) holds with do = 0, R- := cr_ i/2(lnc_ — In do)- 
Since the first addend in (IQ) is now zero, (15.21) shapes into < uv by Remark 15.11 Here uv is the 

maximal acceleration feasible for the robot m, whereas j s the centripetal acceleration required to 

remain on the desired isoline in the “worst case”: the flow is maximal and opposes the robot’s velocity, the 
radius of the isoline is minimal over the range of uncertainty in c and a. 

Let the necessary conditions hold in the enhanced form v > V* and (v + T4) 2 f?Z 1 < vu , where R- := 
min{er_ nc_ — lndo),-R in — (2v + 47 tP*)/m}, and let (15.41) hold with vq := V*. Modulo elementary compu¬ 
tation of the quantities R+ and b n from (15.51) . Proposition 15.21 guarantees that the controller (12.21) succeeds if 
its parameters satisfy (15.511 with 


T/ pa- (v + V t ) 2 

v 0 :=V*, do := 0, A u := u ---, 

vR _ 


u * = u—7T- min Re 2ct - , where u = <py and R+ 

R=R-,R+ 

:= max |cr_ \/2(lnc_ — lnc?o), Rin + (^ v + 47rl4)/u|, 

, |ct 2 — R 2 \ 

b n = max -——. 

R=R-,R + a 2 R 

(T=f7__ ,£7_|_ 


(The last two lines are justified by elementary estimates of min, g [/j| f(z)\ and max, £ [fl_^ + ] TTTzJT'^ 
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6 A Technical Fact Underlying the Proofs of Proposition 14.11 and 
Theorem 14.11 

Lemma 6.1. The following relations hold: 


X(t, r + Tds) = X + ujds + o(ds), 

X(t, r + Nds) = A — Vpds + o(ds) ’ 

(6.7) 

N[t, r + Tds] = N — xTds + o(ds), 

T[t, r + Tds] =T + xNds + o(ds ); 

(6.8) 

N(t, r + Nds) = N + r p Tds + o(ds), 

T(t, r + Nds) = T — r p Nds + o(ds), 

(6.9) 

N[t + dt, r + (dt)] = N — uiTdt + o(dt), 

T[t + dt, r + (dt)] = T + wNdt + o(dt). 

(6.10) 


Proof: Formulas (16.71) are justified by the following: 
A(«,r + Tit) ®> - , v y + %>, = A( t , r) + 


< V ^ T U + D , (E^£i ds <D> A(f , r) _ m Tl d , 


- A 


II vr>„ 

{D"N;T) 
l|V£>|| 
lari 


l|VZA|P 


VD 


, ,, , . , . (VD£ + XD"N;T) 

ds + o(ds) = A (t, r) -- ds 


+ o(ds) = X + cods + o(ds); X(t,r + Nds) = X(t,r) 


(S7 D' t + XD" N ■ N) , to x 

—ds + o(ds) = A — Vpds + o(ds). 


I|V£>|| 


Formulas (16.81) are the Frenet-Serrat equations. Furthermore 

. at! VD(t,r + Nds) AT D"N 
N(> ’ r + Nd ’ } = ||70(t,r + Wd»)|| =N+ M 

{D"N: VD) 


_, , , , , D"N-N(D"N;N) , 

" Vj ° ||V£>||3 + ° ( } “ + [jVDl 

+ o(ds) = N + + o(ds ) JV + TpTds + o(ds). 


This yields the entire (16.91) since N = HjlT.T = —HjlN-, (16.101) follows from the definition of u>. □ 
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